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The static part of the heavy quark potential has been shown to be closely related to the spectrum 
of the rectangular Wilson loop. In particular the lowest lying positive frequency peak encodes 
the late time evolution of the two-body system, characterized by a complex potential. While initial 
studies assumed a perfect separation of early and late time physics, where a simple Lorentian (Breit- 
Wigner) shape suffices to describe the spectral peak, we argue that scale decoupling in general is 
not complete. Thus early time, i.e. non-potential effects, significantly modify the shape of the 
lowest peak. We derive on general grounds an improved peak distribution that reflects this fact. 
Application of the improved fit to non-perturbative lattice QCD spectra now yields a potential that 
is compatible with a transition to a deconfined screening plasma. 



I. INTRODUCTION 

The search for a description of heavy quarkonia at fi- 
nite temperature in terms of a nonrelativistic Schrodinger 
equation with an instantaneous potential has a long his- 
tory. Ever since Matsui and Satz PQ proposed the melting 
of the J/^> particle as signal for the deconfinement tran- 
sition in heavy ion collisions, it has been the goal for the- 
ory to put their phenomenological arguments on a solid 
field theoretical footing. The success of relativistic heavy 
ion experimental groups to indeed measure suppression 
patterns at RHIC and LHC [5] has further invigorated 
interest in the topic. 

Even though initial attempts focussed mainly on model 
potentials [3] , the last decade has seen technical advances 
that allow us to actually derive a potential from the un- 
derlying theory of strong interactions QCD via a system- 
atic coarse graining procedure. Based on the concept 
of effective field theories (EFT) such as pNRQCD g] or 
quantum mechanical path integrals jS] , the static part of 
the in-medium potential can be readily defined from the 
late-time limit of the real-time Wilson loop [6] . 

Using the hard thermal loop approximation, Laine et. 
al. [7] succeeded in calculating the real-time Wilson loop 
to first non-trivial order and obtained a closed expression 
for the potential. They found a real part featuring Debye 
screening and an imaginary part, which was attributed 
to Landau damping [5J. First corrections to this quan- 
tity in an effective field theory approach were derived in 
ref. [5] and the corresponding quarkonium spectral func- 
tions were computed in ref. |10j . 

Questions remain however, such as how the confin- 
ing, i.e. linear part of the potential of the hadronic 
phase behaves when approaching and surpassing the de- 
confinement temperature. Hence we are urged to ex- 
tend the results of perturbative calculations into the non- 
perturbative regime. One path that has proven viable in 
spectral studies of heavy quarkonium [TTJ [T2] is the use 
of lattice regularized QCD, which is amenable to Monte 
Carlo simulations at any temperature. 

As non-perturbative calculations of the Wilson loop 
based on lattice QCD are only possible in Euclidean time, 



we face the challenging task to connect imaginary time 
information with real-time dynamics in the definition of 
the potential. A recent study [T3] suggested that the use 
of Wilson loop spectral functions and the investigation of 
their peak structure allows a non-perturbative extraction 
of both real and imaginary part of the potential. At the 
heart of this procedure is the correct determination of 
the shape of the lowest lying spectral peak. 

From the point of view of the present study, the au- 
thors in [T3] assumed that the scales of early time (bound 
state formation) and late time (potential) physics are 
completely separated. Hence they fitted the spectrum 
with a simple Breit-Wigner and used its peak position 
and width to determine V(r). In the following we wish 
to demonstrate that such an assumption is not valid in 
general and one has to take into account modifications to 
the spectrum that arise from a coupling of the timescales. 
These lead to a skewing and shifting of the Breit-Wigner 
shape. 

After deriving the shape of the lowest lying peak on 
general grounds in section |II A[ we apply the functional 
form to the highest temperature data of ref. [T3] and show 



in section II B that the resulting potential is compatible 
with a scenario of heavy QQ in a deconfined screening 
plasma. The counterintuitive, large rise of the real part 
of the potential found in ref. [13] disappears as a result 
of the improved fit. 



II. HEAVY QUARK POTENTIAL FROM A 
SPECTRAL ANALYSIS 

In order to define a potential for static heavy quarks, 
we rely on the EFT framework. As the constituents of 
the QQ system are heavy [tuq T, Aqcd) and their 
velocities v are small, there exists a hierarchy of scales 
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according to which an EFT description can be con- 
structed. 
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In a first step the hard scale is integrated out to 
give the theory of nonrelativistic QCD (NRQCD), where 
the heavy quarks are described by two-component Pauli 
spinors. One can proceed by integrating out the soft 
scale for the medium degrees of freedom to end up with 
the theory of potential NRQCD (pNRQCD) where color 
singlets and octets are the dynamical fields. 



A. Coupling of scales and spectral shapes 

The matching between QCD and the EFT tells us that 
the time evolution of the Wilson loop [H [5] 



W a (r, t) 



exp 



dx^A, 



in general follows 

id t W a (r,t) = ${r,t)W (r,t) 
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with a time dependent complex function <&(r, t) that 
asymptotes to the singlet potential 



lim $(r,t) = V(r). 

f — >-00 
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To connect this formalism to lattice QCD observables, 
we deploy a spectral decomposition [14] of the Wilson 
loop 



W a (r,t)= / du e** pa(r,u), 



which allows us to extract the potential via 
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We note that if a well defined lowest lying spectral 
peak exists, it will eventually dominate the dynamics in 
the late time limit and thus encodes all necessary infor- 
mation on the potential. In the following we will there- 
fore calculate the shape of the low lying peak, starting 
from the general equations ( |3|4[ ). To this end, we rewrite 
the function $(r, t) — V(r) + <fi(r,t), where the time de- 
pendent part 4>{r, t) vanishes after a characteristic time 
tqQ . Intuitively tqQ is the time needed for the two body 
system of quark and anti-quark to form a bound state 
described by the static potential V(r). It is possible to 
formally solve eq.([3| for t > 0. Bearing in mind that the 
definition of Wrj(r, t) implies Im[V](r) < we obtain 



Wa(r,t) = exp 



i(Re[V}(t)t + Re[a](r, t)) 



\Im[V](r)\t + Im[a](r,t) 



(7) 



Here the function cr(r,t) — f*<fi(r,t)dt is defined as the 
integral over the time dependent part and its asymp- 
totic value <7(xj(r) = c(r,|t| > ^qq) = Jo <fi{r,t)dt. 
From the inversion of eq.(|5| and the positivity condition 
Wn(r, —t) — W n (r, t), we calculate the spectral function: 



/>- ( / . i = ^- J dt exp i (u - Re[V] (r)^t 
iRe[a](r, |t|)sign(t) - |Im[F](r)||t| + Im[a](r, \t\) 



The above expression enables us to write separately the 
actual peak structure coming from the late time physics 
and a background contribution that arises due to the time 
variation of <j)(r, i) 
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The first integral, which runs over the whole time axis can be calculated analytically and will contribute to the well 
defined peak structure encoding the potential. Since we wish to fit the spectra only around the maximum of this 
peak, i.e. where (lj — ReV(r))tqQ <Si 1, we can expand the first term in the second integral exp [i(uj — Re[V A ](r))t] in 
this region. We find that the spectrum can thus be written as 



Pa{r, u) 



r) |Im[V](r)|cos[Re[ CToo ](r)] - (Re[V](r) - w)sm[Re[ CToo ](r)] 



l(r) 



lm[V](r) 2 + (Rc[V](r) - u) 



+c (r) + Cl (r)t QQ (Re[V](r) - u) + c 2 (r)t QQ (Re[V](r) - uj) 2 + 



(9) 



Note that the first term reduces to a simple Breit-Wigner only in the case of dt&(r,t) = 0, which would imply that 
the time- independent potential picture is applicable at all times. In general eq.Q deviates from the simple Lorentian 
through an additional phase and background terms Cj(r), arising from the early time variation of the potential. 
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Even in the region close to the maximum of the peak, 
where all Ci(r) with i > can be ignored, the influ- 
ence of the early time physics modifies the spectral shape 
through c and Re[cr oc ](r) Q In other words, these two 
coefficients have to be considered no matter what fitting 
range is chosen. 

We recover a time independent potential if the spectral 
function Q is plugged back into eq.([6]). Interestingly all 
coefficients Ci{r) as well as contributions from drop 
out in the integratiorj^] To confirm the validity of our 
idea, we also calculated [TS] the spectral function at lead- 
ing order in hard thermal loop resummed perturbation 
theory and observed that it encodes a lowest lying peak 
with exactly the functional form of 

B. Improved fitting on lattice QCD spectra 

The aim of extracting the potential from spectral func- 
tions was to open a window into the non-perturbative 
context of lattice QCD. Monte Carlo simulations allow 
the discrete estimation of correlation functions, such as 
the Wilson loop, in Euclidean time. To connect the real- 
time potential of eq.([6]) and simulated data W\j(r, r, ), we 
have to invert the following Laplace transformation 

/oo 
duj e~ UT Pn (r,uj). (10) 
-oo 

The problem we face is ill-defined, since we wish to ex- 
tract an almost continuous function from a discrete and 
noisy set of points. Fortunately, according to the discus- 
sion following eq.Q, we are interested only in the lowest 
lying peak of the spectrum. This part of the spectrum is 
amenable to analysis via the Maximum Entropy Method 
(MEM) [16- 18 , a form of Bayesian inference. In practice 
its results can depend strongly on the quality of the data, 
i.e. the available signal to noise ratio. Even though the 
position of peaks can be determined reliably, the error in 
the width of spectral structures is more difficult to ascer- 
tain, due to the confluence of systematics and statistics. 

In the following, we will apply the improved fitting 
functions of eq.([9]) to the datasets of [T3]- Spectra are 
available based on quenched lattice QCD data from con- 
figurations with j3 = 6.1, a bare anisotropy of ^ = 
3.2108 and extend 20 3 x 36,24 and 12, which yields 
a„ = 0.097fm and corresponds to the three temperatures 
T = 0.78T C ,1.17T C and 2.33T C respectively [19]. For 
the case of the highest temperature, the number of dat- 
apoints is quite small for the use in MEM, so we provide 



1 This fact was not known in earlier studies e.g. 1131 . which as- 
sumed that a complete separation between early time and late 
time physics holds. 

2 Since eq. (Jsjt is actually a Laurent series, only the principal part 
from the skewed Lorentian, with one negative power, contributes 
to the integrals. 



lattice data as a crosscheck at a more finely spaced set- 
ting of jS = 7, which together with ^ = 3.5 yields a 
lattice spacing of a a = 0.039fm [IT] . 

Fig[l] shows two Wilson loop spectra from the coarser 
lattice at r = 0.19fm and r = 0.39fm, which are fit- 
ted within the shaded region by a naive Lorentian (L), a 
skewed Lorentian (SL, Ci = 0) and the skewed Lorentian 
with up to quadratic background terms (SQL, c L> i = 0). 
We find that the peaks are much better reproduced, once 
the fitting function takes into account at least Co and 
slightly improve with additional c;. The improved fit is 
stable against a change of the fitting range and against 
considering additional Cj. Especially in the right hand 
panel of FigjTJit is obvious that a large background contri- 
bution interferes strongly with a naive Lorentian fit of the 
peak. Note that in the Wilson Loop case at r = 0.39fm 
the signal to background ratio is relatively small, hence 
we only show data up to this distance. 

The values for the real part of the potential Vu { r ) from 
the coarse and fine lattice (the /? — 7 points being shifted 
to account for the different renormalization scale) can 
be found in the left panel of Fig(2j While the naive 
Lorentian fit denoted by (L) shows a very strong rise 
in Re[Vn](r), similar to the one observed in [13] . improv- 
ing the fitting function diminishes this effect significantly. 
We find that at small r, the real part Re[Vb](r) coincides 
within its errorbars, with the temperature independent 
part of the the color singlet free energies and then ap- 
pears to become slightly larger than F 1 (r) at r ~ 0.3fm. 
Most importantly we now understand that the previously 
observed extremely strong rise in the real part can be at- 
tributed to the presence of scale coupling not accounted 
for in the Breit-Wigner fitting function. 

The imaginary part on the other hand does not change 
significantly under the improved fitting. How much of 
the spectral width corresponds to a physical width or is 
induced by an artificial broadening of the MEM is how- 
ever difficult to establish. Without question a width is 
present due to the observed curvature in the Euclidean 
data. Nevertheless until a follow up investigation based 
on e.g. the multilevel algorithm [3U] can simulate the 
Wilson loop at any separation distance with equally small 
relative errors, we suggest to take the right hand side of 
Figj2]as a qualitative result. 

One challenging aspect of the lattice QCD determina- 
tion of the potential is that the Wilson loop in Euclidean 
time becomes more and more suppressed at intermediate 
r, with increasing spatial distance. This suppression in 
turn translates into a decreasing amplitude of the lowest 
lying peak, which hence becomes very hard to extract for 
r > 0.4 just where the non-perturbative range sets in. 

To mitigate these adverse effects, the authors of [T3] 
also considered an alternative observable, i.e. the Wilson 
line correlator WWr, tj) in Coulomb gauge, to define a 
potential Vu (r) in eq. fljj}. The benefit of using this gauge 
fixed quantity is that it shows much less suppression in 
the early r and late r ~ (3 region, which are assumed 
not to contribute to the values of the potential but can 
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FIG. 1. Fits of the quenched (T = 2.337b) lattice QCD Wilson loop spectral function at r = 0.1 9fm (left) and r = 0.39fm 
(right) with a naive Lorentian (L), a skewed Lorentian (SL) and a skewed Lorentian with up to quadratic terms (SQL). Fitting 
with a naive Lorentian is inadequate already at the second smallest separation distance. We find that skewing alone does not 
remedy the fit, while including up to quadratic terms allows us to reconstruct the spectral peak very accurately. Note that the 
naive Lorentian fit systematically overestimates the value of the real-part of the potential for r > O.lfm. 
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FIG. 2. (left) Real part of the lattice QCD heavy quark potential Vfo(r) at T = 2.337b. For comparison we show the 
color singlet free energies at T = 2.337b and T = 0.787b (gray solid, dashed lines). We find that the extraction of the real 
part changes significantly if improved fitting functions are deployed. The strong rise of Re[Vb] vanishes, so that its values 
become only slightly larger than the color-singlet free energies. For completeness we show (right) the reconstructed values of 
the imaginary part of the potential from the fit of the spectral width. 



complicate their extraction [BJ. 

The shape of the spectra from the Wilson line correla- 
tor for the coarser lattice tell us that the background con- 
tributions are significantly reduced compared to the Wil- 
son loop case. We are thus able to identify a well defined 
peak up to distances of r ~ lfm. At small separation dis- 
tances r < 0.4fm, the naive fitting with a Lorentian and 
the improved fitting functions yield a very similar peak 
position. Only by going to larger separation distances 
the improved fitting gives significantly better results. 

Turning to the reconstructed values of the potential 
V[|(r) from the coarse and fine lattice as shown in Figj3j 
we find that the real part shows a behavior much closer 
to intuition as was found without improved fitting func- 
tions. Instead of a perpetual linear rise, Re[Vj|](r) grows 
only until distances of r ~ 0.45fm, where it begins to 
flatten off to a value slightly larger than the color-singlet 



free energies. Note that the real part of Vj| (r) and Vrj(r) 
mostly agree within their error bars in the range where 
Vrj(r) can be determined. The imaginary part Im[Vj|](r) 
is however smaller than Im[Vbj](r), which could be a sign 
of better acuity in the spectral peaks due to a higher 
signal to noise ratio. 



III. CONCLUSIONS AND OUTLOOK 

We have shown that the coupling of early and late 
time scales leads to spectral structures different from a 
Breit-Wigner, which have to be taken into account to ex- 
tract the heavy quark potential reliably. After deriving 
improved fitting functions in ^ , we applied them to the 
non-perturbative Wilson loop spectra and found that the 
fits reconstruct all spectral functions excellently, which is 
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FIG. 3. (left) Real part of the lattice QCD heavy quark potential Vji(r) at T = 2.33Tc. For comparison we show the color 
singlet free energies at T — 2.33Tc and T = 0.787c (gray solid, dashed lines). We find that the extraction of the real part 
between different fitting functions only starts to differ for values above r > 0.4fm, due to the smaller background contributions 
in p\\(r, uj). The real part Re[Vji](r) grows to values slightly larger than the color singlet free energies before it asymptotes to 
a constant value at r > 0.5fm. (right) Values of the imaginary part of the potential from a fit of the spectral width. 



by far not the case for the Lorentian. With our improved 
fitting, we are able to determine a real- and imaginary 
part of the potential, which is compatible with the pres- 
ence of a deconfined and screening quark-gluon plasma. 

The error bars for the extracted values of the potential, 
especially in the case of the Wilson loop, are still rela- 
tively large. Thus to make more quantitative statements 
in the future, efforts need to be continued to increase the 



signal to noise ratio of the Euclidean correlator data in 
order for the MEM to improve. 
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